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ABSTRACT 

We present new results from BRAVA, a large-scale radial velocity survey of the Galactic bulge, using 
M giant stars selected from the Two Micron All Sky Survey catalog as targets for the Cerro Tololo 
Inter-American Observatory 4m Hydra multi-object spectrograph. The purpose of this survey is to 
construct a new generation of self-consistent bar models that conform to these observations. We report 
the dynamics for fields at the edge of the Galactic bulge at latitudes b = — 8° and compare to the 
dynamics at b = —4°. We find that the rotation curve V(r) is the same at b = —8° as at b = —4°. That 
is, the Galactic boxy bulge rotates cylindrically, as do boxy bulges of other galaxies. The summed line 
of sight velocity distribution at b = —8° is Gaussian, and the binned longitude-velocity plot shows no 
evidence for either a (disk) population with cold dynamics or for a (classical bulge) population with 
hot dynamics. The observed kinematics are well modeled by an edge-on TV-body bar, in agreement 
with published structural evidence. Our kinematic observations indicate that the Galactic bulge is a 
prototypical product of secular evolution in galaxy disks, in contrast with stellar population results 
that are most easily understood if major mergers were the dominant formation process. 
Subject headings: Galaxy: bulge - Galaxy: kinematics and dynamics - stars: late-type - techniques: 
radial velocities 



1. INTRODUCTION 

The central bulge of the Milky Way is our nearest ex- 
ample of a spheroidal population, with M31 100 times as 
distant. The Galactic bulge stellar population by far is 
the nearest bulge population that can be studied in star- 

1 Visiting Astronomers, Cerro Tololo Inter- American Observa- 
tory (CTIO). CTIO is operated by AURA, Inc. under contract to 
the National Science Foundation. 



by-star detail. Being close enough to permit the study 
of radial velocity, proper motion, and composition for in- 
dividual stars, as well as turnoff age for the population, 
the bulge/bar population offers us an unprecedented op- 
portunity to test dynamical and formation models for 
bulge systems. This is a unique perspective unavailable 
in the study of extragalactic bulges, which can only be ex- 
amined from their integrated light, and may change the 
way we think about the formation of these structures. 
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The Bulge Radial Velocity Assay (BRAVA) exploits this 
unique opportunity with a large scale radial velocity sur- 
vey of the Galactic bulge/bar population. 

The boxy morphology of the central bulge is easily seen 
in the image of the Galaxy produced by the COBE satel- 
lite (Weiland et al.1994; Dwek et al. 1995; Arendt et al. 
1998) and subsequent models have solidified the interpre- 
tation that we view an edge-on bar (e.g Zhao et al. 1996; 
Hairier et al. 2000). Exploiting our proximity, we know 
that the bulge/bar is old (>10 Gyr) from the modeling 
of its main- sequence turnoff (Ortolani et al. 1995; Kui- 
jken & Rich 2002; Zoccali et al. 2003; Picaud & Robin 
2004; Clarkson et al. 2008). Observations of individual 
stars at high resolution yields evidence of a enhance- 
ment (McWilliam & Rich 1994; Fulbright, McWilliam, 
& Rich 2007; Lecureur et al. 2007) that is modeled to 
imply a short formation timescale of < 1 Gyr (Ballero et 
al. 2007). And while there are hints of an intermediate 
age population in the bulge at lower Galactic latitudes 
(van Loon et al. 2003) the bulk of the bulge lacks a con- 
vincingly demonstrated intermediate age evolved stellar 
component (Frogel & Whitford 1987). These studies re- 
inforce the widely accepted paradigm that the bulge pop- 
ulation was formed both rapidly and early. 

Combes & Sanders (1981) were the first to suggest 
that galaxy bars heat themselves in the vertical direction 
and look boxy when seen edge-on. Additional observa- 
tions and modeling have confirmed the general picture 
that boxy and peanut bulges are not spheroidal merger 
remnants (hereafter "classical bulges") but are in fact 
edge-on bars (e. g., Combes et al. 1990; Athanassoula 
& Misiriotis 2002; Athanassoula 2005; see Kormendy & 
Kennicutt 2004 for a review). However, in our Galaxy, 
it is difficult to reconcile the early and rapid star for- 
mation implied by measurements of bulge star ages and 
a-element overabundances with the protracted star for- 
mation that might be expected for a secularly growing 
bar (Kormendy & Kennicutt 2004, see Section 8.1). Free- 
man (2008) emphasizes that "in the bar-buckling insta- 
bility scenario for generating boxy bulges, ... the bulge 
structure may be younger than its stars, which would 
originally have been part of the inner disk." So the old 
age and rapid star formation history implied by observa- 
tions of the Galaxy's boxy bulge may not be a problem 
for our interpretation that this boxy structure was built 
out of the disk. 

Further tests of whether the Galaxy's boxy bulge really 
is a edge-on bar would therefore be valuable. Classical 
bulges are observed to rotate more slowly at increasing 
height above the disk plane, whereas boxy bulges in other 
galaxies rotate cylindrically, with V(r) essentially inde- 
pendent of height above the disk plane (e. g., Kormendy 
& Illingworth 1982; Jarvis 1990; Shaw 1993; Falcon- 
Barroso et al. 2004; see Kormendy & Kennicutt 2004 
for a review). Supporting the above interpretation, TV- 
body models of bars also rotate cylindrically when viewed 
edge-on (Combes et al. 1990; Fux 1997, 1999; Zhao et 
al. 1996; Athanassoula & Misiriotis 2002; Athanassoula 
2005). To test our interpretation, we ask: does the boxy 
bulge of our Galaxy rotate cylindrically or not? 

There is also presently some debate as to whether the 
bar continues to dominate the population at 1 kpc from 
the bulge. Zoccali et al. (2008) suggest that the emer- 
gence of a metal poor classical bulge population at higher 
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Fig. 1. — Observed BRAVA fields up to 2007 August overplotted 
on the COBE 2 fim image (Launhardt et al. 2002). Green circles 
represent fields observed in 2005, red circles are fields observed 
in 2006, and blue circles are fields observed in 2007A (b = —4°) 
and 2007B (b = —8°). The yellow circle is our calibration field, 
observed with the same fiber configuration every year. Circle sizes 
correspond to the 40' field of view of the instrument. The black 
square is a field observed over multiple years. For a distance of 8 
kpc, 1° corresponds to 140 pc. 

Galactic latitudes might explain why an abundance gra- 
dient is observed in the bulge, as a secularly evolved bar 
should not have an abundance gradient. Further, the bar 
has been solidly demonstrated at b = —4° (e.g. Stanek et 
al. 1997) but is not convincingly detected more distant 
from the plane. A large sample radial velocity survey 
at b — —8° offers a sensitive test of whether a classical 
bulge population is present. 

Until recently, we have lacked a kinematic survey of 
the bulge large enough to seriously test these bar and 
bulge dynamical models, and thus influence ideas about 
bulge evolution. The kinematic studies of bulge red giant 
branch (RGB) stars from BRAVA (Howard et al. 2008) 
at b = —4° and bulge red-clump giants (Rangwala et al. 
2009), show remarkable agreement in the dispersion and 
rotation curve. Measurements of the kinematics of the 
b = — 8° strip are especially important, since it gives us 
an opportunity to investigate whether the bulge rotates 
cylindrically. Here we report results from a continuing 
survey based on red giants, BRAVA, which comprise the 
bulk of the 2.4 /im light of the bulge. Early results from 
this ongoing survey are given in Rich et al. (2007b) and 
Howard et al. (2008), which report on the 2005-2007A 
observations. This Letter reports on the 2007B observa- 
tions of the b = —8° major-axis strip. 

2. OBSERVATIONS AND SPECTROSCOPY 

We report on 12 new bulge fields observed in 2007B, 
sampling the major axis of the bulge in two degree in- 
tervals at b = —8° (Figured]). We obtained radial veloc- 
ities for ~100 red giants per field. We used the Hydra 
multi-object spectograph at Cerro Tololo Inter- American 
Observatory (CTIO); the observational setup and target 
selection were as described in Howard et al. (2008). Be- 
cause sources that satisfy our selection criteria are rare 
at b = —8°, we extended our selection limit by ~0.5 mag 
in K compared to our selection region for the 2005-2007A 
data (see Figure 9 of Howard et al. 2008). Although we 
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are observing fainter stars in 2007B than in 2005-2007A, 
our selection region is still brighter than the "red-clump" 
region and samples the same bulge RGB population. The 
2007B radial velocity standards (HD 203638, HD 207076, 
and HD 218541) return individual stellar velocities that 
agree to better than 2 km s _1 , on average. A final ve- 
locity was obtained by taking the error- weighted average 
of the separate velocity measurements. One field was 
observed in all three years (2005, 2006, and 2007B) at 
(Z, b) = (6°, —4°), with the same fiber configuration (i.e., 
the same stars in that field were observed each year). 
Sky conditions during the 2007B night when this field 
was observed gave less-than-ideal results that are not in- 
dicative of the 2007B data in general. Velocity measure- 
ments show a larger rms scatter in velocity differences 
(—7 km s _1 ) as compared to 2005/2006 data, with an 
average offset of —2 km s _1 for all stars in that field. 
Since the offset is less than the rms scatter, we consider 
the 2005/2006/2007B data sets to be in good agreement 
and adopt —7 km s _1 as our error estimate for individual 
stellar velocities for all of the data obtained in 2007B. 

2.1. Color/Magnitude bias 

Observations of the bulge at b = — 8° have the po- 
tential to probe the bulge/halo boundary. We examine 
these b = — 8° separately from the higher latitude fields 
presented in Howard et al. (2008). To investigate the 
possibility of color and/or magnitude bias in our sam- 
ple, the bulge RGB stars from the b = — 8° fields are 
summed and yield an apparent Gaussian distribution, 
with < V G c >= -9-1 ± 2.7 km s -1 and a = 94.4 ± 1.9 
km s _1 (Figure [2]). The longitude- velocity (l-v) plots for 
both the b = -4° (see Figure 20 of Howard et al. 2008) 
and b = — 8° (Figure [3]) major- axis strips show roughly 
trapezoidal "envelopes" that show no hints of contam- 
inating hot (spheroid), or cold (disk), components. In 
order to determine if there are color /magnitude biases 
in our sample, we divide the sample by color and mag- 
nitude (see Howard et al 2008 for details) and employ 
the two-sided Kolmogorov-Smirnov (KS) test. As with 
the 2005-2007A data (Howard et al. 2008), we can state 
that the 2007B color /magnitude-segregated populations 
are drawn from the same distribution, rejecting the null 
hypothesis with 96% confidence. At b — —8°, we are at 
the point where the bulge density begins to drop sharply; 
our aim is to search for subpopulations in these regions, 
and the signature of the inner halo (C. D. Howard et al. 
2009, in preparation). 

3. EVIDENCE OF CYLINDRICAL ROTATION 

We now discuss the Galactic boxy bulge's kinematics 
at —560 and —1120 pc below the disk plane. Figured 
shows the dispersion profile and rotation curve for the 
major axis strips at b = —4° and —8°. Despite the more 
sparse sampling at b = —8°, it is clear that the rotation 
curve here is indistinguishable from that at b = —4°. We 
therefore confirm that the Galaxy's bulge rotates cylin- 
drically. This is further strong evidence that it is an 
edge-on bar. It is difficult to conclusively argue whether 
the rotation curve flattens like that observed at b = —4°, 
however, the inclusion of data obtained at this latitude 
in 2008 will be able to determine if this is the case (C. 
D. Howard et al. 2009, in preparation). 
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Fig. 2. — Histogram of all bulge RGB star velocities from the 
fields at b = —8° (~1200 stars), in galactocentric velocity (Vgc)- 
As with the 2005-2007A data, the co-added sample is consistent 
with a single kinematic population that is normally distributed 
and that has negligible skew and kurtosis. Bin size is 25 km s _1 . 
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Fig. 3. — Longitude-velocity (l-v) plot for the 2007B bulge sam- 
ple at b = —8°; the data are binned to 1° in longitude and 10 
km s — 1 in galactocentric velocity. There is no evidence for either 
a cold, rapidly rotating disk, or for a dynamically hot, more slowly 
rotating (more classical bulge-like) population. The velocity dis- 
persion drops from that at —4° and is also lower than the halo. 
We conclude that the dynamical bar population is predominant at 



3.1. Model Comparison 

A useful comparison of model and data involves the 
disentangling of the different stellar components of the 
model, i.e., spheroid versus disk/bar components, and 
comparing their kinematics to our observations. We com- 
pare our data to an TV-body bar model (Fux 1999). The 
model is constructed from a composite three-dimensional 
symmetry-free TV-body and hydrodynamics code which 
follows the constituent particles at a higher resolution 
than previous models and includes a gas component 
which reproduces the CO and HI distributions in the 
(l,b,V) space (Fux, 1999). It includes -1.5 x 10 5 parti- 
cles of gas, —1.3 x 10 6 particles in the stellar disk/bar, 
—6 x 10 5 particles in a nucleus-spheroid (representing 
a spheroidal nucleus and a stellar spheroid outside the 
bar region), and — 1.7 x 10 6 of dark halo particles. In 
Figure [51 we compare the BRAVA data to the disk/bar 
and nucleus-spheroid components of the model named 
'cl0t2066' described in Fux (1999). 
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Fig. 4. — Velocity dispersion profile (top) and rotation curve 
(bottom) for the b = —4° and —8° strips. The velocity dispersion 
drops as expected, while the rotation curve is the same for both 
latitudes, consistent with cylindrical rotation. The b = —8° data 
are also consistent with solid-body rotation, and at first inspection 
does not appear to show the flattening observed at b = —4° . More 
observations are needed to confirm this finding. 

As one moves further in latitude from the Galactic disk 
plane, it is evident that the nucleus-spheroid component 
of the model is hotter and contributes a large velocity dis- 
persion. Despite this, the BRAVA data show remarkable 
agreement with the disk/bar component of the model, 
with a relatively flat dispersion profile at ~70 km s _1 
contrasting with the spheroid dispersion of ~120 km s _1 . 
Despite the poorer sampling of fields in 2007B, the data 
show good agreement with the rotation curve as well, and 
suggest that the turnover seen in the rotation curve at 
b = —4° is not evident at b = —8°. Of course, this can- 
not be confirmed until the rest of the b = —8° major-axis 
strip data are included. 

4. CONCLUSIONS 

Following the first successful models of the boxy mor- 
phology of the 2.4 /im light (Blitz & Spergel 1991) the 
case for the deprojected bar has become more convinc- 
ing (e.g., Bissantz & Gerhard 2002 models of the COBE 
light) and is reinforced by distances of red-clump stars 
(e.g., Stanek et al. 1997; Babusiaux & Gilmore 2005). 
The boxy outer isophotes and our measurement of cylin- 
drical rotation at b = — 8° are also consistent with a bar 
viewed edge-on; this is one proposed class of pseudobulge 
as defined by Kormendy & Kennicutt (2004). We empha- 
size that the rotation curve is strikingly inconsistent with 
a slowly rotating, dynamically hot, classical bulge popu- 
lation. The excellent agreement of the disk/bar compo- 
nent of the Fux (1999) TV-body model with the BRAVA 
data further supports our interpretation of the Galac- 
tic pseudobulge as an edge-on bar. We assert that the 
Galactic bulge formation scenario does not adhere to the 
hierarchical clustering paradigm (which views the bulge 
as a merger remnant) but rather represents the secular 
evolution of isolated disks. 

A range of TV-body models suggest that the timescale 
to grow a bar from a massive disk is very short, as few 
as 2-4 orbital times (Combes et al. 1990; Raha et al. 
1991; Sellwood & Wilkinson 1993, Fux 1997, 1999). In 
our Galaxy, this corresponds to ~ 10 8 yr at 1 kpc ra- 
dius. Note that other acceleration mechanisms (e.g. ver- 
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Fig. 5. — Upper: results for the b = —4° major-axis strip includ- 
ing dispersion (top) and rotation field (bottom), compared with 
Fux (1999) model components. Lower: the b = —8° major-axis 
strip dispersions (top) and mean field velocity (bottom), compared 
with Fux (1999) model components. The difference in velocity dis- 
persion between the spheroidal and disk/bar components is strik- 
ing, and shows that BRAVA data are well modeled by an TV-body 
bar. The rotation speed does not drop with increasing distance 
from the plane (cylindrical rotation); such cylindrical rotation is a 
signature of a boxy bulge. 

tical resonances; Pfenniger & Norman 1990) appear to 
work on longer timescales — at least ~ 10 orbital times 
to achieve vertical thickening. The observed a enhance- 
ments measured for pseudobulge giants (McWilliam & 
Rich 1994; Rich & Origlia 2005, 2007a; Zoccali et al. 
2007; Fulbright et al. 2007) would appear to require 
massive star supernovae, with enrichment proceeding on 
10 8 yr timescales (Ballero et al. 2007). One may cau- 
tion that much larger samples of composition measure- 
ments, across the entire pseudobulge, are required before 
the chemical composition becomes a very strong con- 
straint on formation scenarios. However, present sam- 
ples of bulge giants are now large enough to permit us to 
conclude that a enhancement is a characteristic of pseu- 
dobulge stars. 

4.1. Discussion 

If the Galactic "bulge" is a pseudobulge that secularly 
evolved via dynamical processes, it is difficult to under- 
stand how a metallicity gradient, as reported by Zoccali 
et al. (2008), could arise. Proctor et al. (2000) find 
a strong abundance gradient in NGC 4565, which is an 
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edge-on spiral with a boxy bulge. The Milky Way falls 
near NGC 4565 on the (V max /(j) versus e plot (Howard 
et al. 2008, Rich et al. 2007b), and the boxy bulge of 
NGC 4565 also shows cylindrical rotation (Kormendy & 
Illingworth 1982). An abundance gradient in our boxy 
pseudobulge might arise if there were an old metal-poor 
spheroid component spatially coincident with the bar. 
We argue that Figures 2 and 3 all but rule out a dy- 
namically hot subcomponent that one might assign to 
such a metal-poor spheroid or classical bulge. A minor 
merger might also leave an abundance gradient but would 
be expected to have left a dynamically cold component 
behind or possibly large population of evolved interme- 
diate age stars (e.g., carbon stars or luminous M giants 
as are seen in the LMC bar) if the ingested system had 
an age range. No such stellar population is known in the 
boxy pseudobulge at high Galactic latitudes, and large 
numbers of carbon stars, for example, would have been 
easily detected. An abundance gradient might also arise 
if the bar were imprinted rapidly during the violent star- 
forming phase of the bulge; the relative confinement of 
metal-rich stars to the plane would argue for some dissi- 
pative process accompanying enrichment. Soto, Rich, & 
Kuijken (2007) find that the bar is most strongly evident 
in the proper-motion dynamics of the metal-rich popu- 
lation, although the kinematics of stars in our b = — 8° 
fields suggest that virtually all stars in our sample belong 
to the bar. 

We remain puzzled as to how the boxy pseudobulge 
could have evolved on a rapid dynamical timescale, yet 
also have an abundance gradient. If the stars are formed 
and are in a massive disk once the buckling process starts, 
it is difficult to understand how any abundance gradient 
can be imprinted on the resulting bulge, since the TV-body 
models accelerate points undistinguished by any physi- 
cal property. We believe that our data are now strong 
enough that one has great difficulty proposing a boxy 



pseudobulge that transitions into a metal-poor classical 
bulge at b ~ —8°: the required dynamically hot, slowly 
rotating, population is not observed. While such a tran- 
sition would help to explain the gradient, our BRAVA se- 
lection criteria at b = —8° should include such stars; the 
dynamics of the population are remarkably uniform and 
leave little room for this hypothetical population. We 
must conclude, then, that the old population at b = —8° 
is dominated by the same box/peanut pseudobulge pop- 
ulation observed at Galactic latitudes closer to the disk 
plane. Even if two types of bulges (both "classical" and 
"box/peanut") coexist in these low latitude fields, it re- 
mains problematic as to how there can exist an abun- 
dance gradient with such uniform kinematics, but dissi- 
pation during the chemical enrichment process offers an 
interesting route toward a solution. The BRAVA survey, 
however, is finding no indication that a classical bulge 
population sets in 1 kpc below the plane. 
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